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SUMMARY 


An investigation of the effect of leading-edge thickness on the flow 
over flat plates with square and cylindrical blunting was conducted at a 
Mach number of 4 and free- stream Reynolds numbers per inch of 2380 and 
6600. Surface pressures were measured on a series of models whose leading- 
edge thicknesses ranged from 0.25 to 1 inch. Heat-transfer rates were 
measured from a flat plate which was blunted by a 1- inch-diameter cylin- 
drical leading edge. All tests were performed with the instrumented 
surfaces at zero angle of sweep and zero angle of attack. 

For the test conditions, the bow shock wave was detached and leading- 
edge shape had no effect on surface pressures aft of two leading-edge 
thicknesses. The surface pressures could be predicted by a combination 
of shock-wave boundary- layer interaction theory and blast wave theory. 

This combination applied equally well to similar data of other investiga- 
tions. An empirical expression for local Reynolds number at the boundary- 
layer edge was found to correlate both the present data and data from 
other investigations covering a wide range of conditions. The local 
Reynolds number per inch was found to be lower than free-stream Reynolds 
number per inch, nearly constant for the test length, and to have negli- 
gible dependence on leading-edge bluntness. This reduction depends on 
the square root of the ratio of total pressures across the normal bow 
shock wave. 

The local Nusselt number was found to depend only on the local 
Reynolds number for the present tests, and is predicted by the familiar 
Pohlhausen flat-plate theory. As compared to the sharp condition, blunt- 
ing the leading edge of flat plates, with consequent reduction of local 
total pressure, was found to increase the heat-transfer coefficients in 
the region where surface static pressures were high and to reduce the 
coefficients where the surface static pressures approached the free-stream 
value . 
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INTRODUCTION 


The presence of the boundary layer on a body effectively changes body 
shape. The effect on a slender body may be sufficiently large to produce 
a measurable strengthening or even detachment of the bow shock wave at 
supersonic speeds. Also, as speed is increased, interactions occur between 
the boundary layer and the shock wave, because the shock wave "wraps" 
closely back over the body. These two phenomena, shock strengthening 
and "wrapping" effect, change the flow field from that occurring in the 
usual low supersonic speed, high Reynolds number laminar flow over a 
plate . 

A theoretical study of the hypersonic viscous flow over a sharp- 
nosed plate (ref. 1) brought out a parameter to describe phenomena occur- 
ring in the region of "weak" interaction between the boundary layer and 
the leading-edge shock wave. This interaction parameter is a special 
combination of the usual flow parameters, Reynolds number and Mach num- 
ber. The region of strong interaction, near the leading edge, was inves- 
tigated theoretically in reference 2 where this same interaction parameter 
was used. Correlation of surface static pressures measured on flat plates 
with square leading edges in hypersonic flow was obtained in reference 3 
using this interaction parameter, but the magnitude of the pressures was 
greater than that predicted by the theories of references 1 and 2. The 
discrepancy between the theory and experiment was ascribed in reference 3 
to leading-edge bluntness effects. Also, high static pressures measured 
by Kendall (ref. 4 ) on a flat plate with a very sharp leading edge were 
correlated by the interaction parameter but not satisfactorily predicted 
by the interaction theories of references 1 and 2. 

In reference Hammitt, Vas, and Bogdonoff formed an empirical 
relation to explain the leading-edge thickness contribution to the surface 
pressures measured on a series of blunted flat plates. A linear combina- 
tion of the viscous interaction parameter and an empirical "inviscid" 
parameter was postulated. They obtained fair correlation of their own 
tests using this empirical relation; however, they were unable to predict 
the results of the experiments described by Bertram in reference 6. 

More recently, blast wave theory has been applied to the inviscid 
problem by Lees (refs. 7 and 8) to obtain a functional relation for the 
inviscid pressure term. Cheng and Pallone (ref. 9 ) , following similar 
lines of reasoning, obtained an inviscid pressure term similar to that 
of Lees. Both differ from the empirical relation of reference 5. 

Lees, reference 10, analyzed the influence of the leading-edge shock 
wave (from blunt and sharp leading edges) on the laminar boundary layer 
at hypersonic speeds. For sharp leading edges he showed the effects of 
a highly curved attached leading- edge shock to be carried back over a 
large portion of the surface and to alter appreciably the temperature and 
vorticity at the outer edge of the boundary layer. By' considering detached 
shocks for a blunt leading edge, Lees predicted the change of induced 
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pressure with interaction parameter as measured ty Hammitt and Bogdonoff 
and described in reference 3 . However, the magnitude of the pressure 
was not predicted. 

The purpose of the research described in this paper was to study 
the aerodynamic and heat-transfer processes occurring over the surface 
aft of the leading edge of flat plates with detached shock waves. It 
was hoped that the results obtained with the low-density test conditions, 
which provide an aerodynamic magnification of the leading-edge region, 
would clarify the effect of size and shape of leading edge on the 
fundamental processes occurring there. 

Experimental pressure-distribution and heat-transfer data were 
obtained from blunted flat plates in a rarefied gas stream. The results 
are compared with the results of various other similar investigations. 
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SYMBOLS 


a 

b 


B.L. 


c 


C w 

c i 

C D 

d 

h 

I 

Io 

k 

k f 

K 


exponent defined in equation ( 8 ) 
constant defined in equation (4b) 
boundary layer 

factor defined in equation ( 5 b) 
specific heat at constant pressure, Btu/lb, °F 
constant in linear relation between viscosity and temperature 
constant defined in equation ( 5 b) 
drag coefficient based on frontal area 
diameter or thickness, ft 
heat-transfer coefficient, Btu/ft 2 , hr, °F 
blast wave pressure term. 


Hx , 2 


(x/d) 


2/3 


modified Bessel function, zero order, first kind 
thermal conductivity, Btu/hr, ft 2 , °F/ft 
film thermal conductivity, Btu/hr, ft 2 , °F/ft 
constant defined by equation ( 8 ) 
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Ko 

modified Bessel function, zero order, second kind 

M 

Mach number, dimensionless 

N 

number of molecules per unit volume, l/cu ft 

Nu x 

Vjy- 

Nusselt number, — , dimensionless 
k 

P 

pressure, lb/sq ft 

Pr 

3600 pc- 

Prandtl number, — , dimensionless 

k 

q. 

heat flow, Btu/hr 

r B 

outer radius of annular space around heat-transfer element, ft 

r m 

radius at minimum temperature point on film, ft 

r 

radius of point on film, ft 

Re <=od 

P ii d 

Reynolds number, , dimensionless 

Poo 

Re x 

Reynolds number , dimensionless 

r" 

s 

area defined by rp, sq ft 

tf 

film thickness, ft 

T 

temperature, °R 

To 

temperature of test-section surfaces, °R 

u 

velocity in x direction, ft/sec 

v m 

most probable molecular speed, ft /sec 

X 

coordinate length, ft (see fig. 3) 

y 

coordinate length, ft (see fig. 3) 

a, 3 

empirical constants 

7 

ratio of specific heats, dimensionless 

5 

boundary- layer thickness, ft 

5 

parameter defined in equation (B9) 
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e emissivity, dimensionless 

k Boltzmann constant, 7.28xl0" 27 , Btu/molecule, °R 

p viscosity, lb/sec ft 

p density, lb/cu ft 

a Stefan-Boltzmann constant, 1.73x10" 9 , Btu/ft 2 hr, °R 4 

X interaction parameter, 

w exponent of viscosity-temperature law defined in equation (B4) 

Subscripts 

b blunt 

d characteristic dimension for quantity involved 

p heat-transfer element quantity 

s sharp 

ti total quantity in free stream ahead of all disturbances 
t 2 total quantity behind leading-edge normal shock wave 

t 3 total quantity behind probe normal shock wave 

tg local total quantity 

w body surface quantity 

x characteristic length for quantity involved 

T) "effective" recovery condition 

oo quantity evaluated at free-stream condition ahead of all disturb- 

ances 

Superscript 

( ) 1 T 1 method 
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DESCRIPTION OF EQUIPMENT AND TEST METHOD 
Wind Tunnel 


The tests were conducted in the Ames 8-inch low-density wind tunnel. 
This wind tunnel is an open- jet nonreturn type tunnel. Air was used as 
the test gas. The 8-inch tunnel is a scaled up version of the low-density 
wind tunnel described in reference 11. A five-stage set of steam ejectors 
is used to produce the main flow. The axisymmetric nozzle was designed 
by the method described in reference 12. The nozzle was constructed of 
shim stock of varying thickness and alternate shims were removed to per- 
mit boundary- layer removal as described in reference 13 . The design Mach 
number was 4 through the stream-static pressure range of 100 to 300 microns 
of mercury absolute. The boundary layer is removed by steam ejectors 
operating in parallel with the main drive set. The physical arrangement 
of the nozzle and test section is shown in figure 1. 

Impact pressure surveys indicated that no strong shock waves were 
present in the nozzle and jet when the expansion ratio across it was 
properly set and controlled. Additional surveys were made in a plane 
normal to the stream direction 1-1/4 inches downstream of the nozzle exit. 
The static pressure of the stream was obtained by measuring the nozzle 
wall pressure at a point 2 inches upstream of the exit plane of the noz- 
zle. This method of obtaining stream static pressure has been described 
in reference 12. 

A typical Mach number distribution obtained from these measured 
quantities is shown in figure 2. The Mach number was calculated in two 
ways, (l) from measured impact pressure and static (wall tap) pressure 
together with Rayleigh’s pitot fo'rmula (circular symbols), and (2) from 
measured impact pressure and upstream reservoir pressure (total head) 
using the assumptions that the flow through the nozzle was isentropic and 
that the probe produced a normal shock wave (square symbols). Good agree- 
ment was obtained between the two methods of obtaining Mach number over 
the range of pressure levels used in the investigation. Therefore, the 
assumption that the flow through the nozzle was isentropic appears to be 
reasonable . 

Table I presents the actual usable stream diameter and Mach number 
for two test-section static pressures. 


Pressure Models 


Four pressure models were tested. Three two-dimensional wedges 
of 5° total included angle were instrumented for measurement of surface 
pressures. The leading edges were cylindrical of 1, l/2, and l/4 inch 
diameter. A 1- inch- thick slab with a square leading edge was similarly 
instrumented. The models were sufficiently wide (6 in.) to span the 
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usable stream completely. The length of afterbody from shoulder point 
(tangent point) was 5 inches. Pressure orifices, l/32-inch diameter, 
were installed in the upper surface at intervals of l/l6 to 1/2 inch 
along the midspan plane. The spacing was smallest near the leading edge. 
An orifice was installed at the nose of each body on the axis of symmetry. 

The surface orifices were connected to a pressure manifold by stain- 
less steel tubes. A pressure switch was used to select and connect any 
desired orifice to the pressure-sensing device. An oil manometer was 
used to measure the orifice pressures. The manometer fluid was diffusion 
pump oil which had been well outgassed. 

Surveys were made in the flow field adjacent to the body with an 
impact tube. The tube end was flattened to an over-all height of 
about O.O3O inch to minimize disturbances in the y direction. The 
pressures from this probe were likewise measured with an oil-filled 
U-tube manometer. 


Heat-Transfer Model 


The heat-transfer model, figure 3, was similar in shape to the 
cylindrical nosed pressure distribution models described above. The 
model was constructed of copper with a leading-edge diameter of 1 inch. 

The model was 5 inches in length from the shoulder to the base. It was 
6 inches wide. Cylindrical heaters installed in the body, one on each 
side, were used to maintain the desired model temperature. 

At various locations along the model length, heat-transfer elements, 
hereafter called heat plugs, were installed as seen in figure 3. The 
heat plug is a copper spool, l/8-inch diameter and l/8-inch long, inserted 
in a hole in the model shell with a 1/64- inch air gap between the plug 
and the body. The spool is supported by a Bakelite stem. The surface 
of the plug was machined to the contour of the plate. An electrical 
heating coil was wound on the plug to supply heat to the plug. A differ- 
ential thermocouple indicated the temperature difference between the plug 
and the model body. Thermocouples were mounted in the body near the heat 
plug to indicate temperature of the body. 

A 0.003-inch-thick plastic film was placed over the surface of the 
model to seal the gap between the plug and the body from the air stream. 
The air gap was then vented to the hollow portion inside the model and 
thence to constant static pressure of the test chamber. Thus, the heat 
loss due to conduction through the air gap was minimized because of the 
presence of a quiescent layer of low-pressure air around the test plug. 
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Test Method 


The heat-transfer rate from the surface of the plug to the air stream 
was determined as a function of the plug temperature. A test point was 
obtained by heating the model and the test plug to the same uniform tem- 
perature, and measuring the plug heater current for this steady-state 
condition. 

A series of tests were made, with no air flow through the tunnel, 
at pressures of approximately 0.1 micron of mercury absolute and of 100 
to 300 microns, to obtain the radiation and conduction heat losses. The 
variation of the sum of these losses with pressure was within the scatter 
of the data. This total loss was then treated as a tare to be subtracted 
from the gross heat input to the plug obtained in each test with air flow. 
The magnitude of the tare was approximately 10 percent of the total heat 
input at the highest rate and about 60 percent of the total heat input 
at the lowest rate. The lowest heat rate occurred near the rear of the 
plate at the lowest stream pressure. 

At a given flow condition, net heat input to the plug was obtained 
at a series of plug temperature levels ranging 20 ° to 50 ° P above stagna- 
tion temperature. The net heat flow to the plug, which is the heat flow 
to the stream, was divided by the heat-transfer area and plotted as a 
function of the difference between the plug temperature and the stagnation 
temperature of the stream. The slope of this curve is proportional to 
the heat-transfer coefficient. The heat-transfer area was taken as 
slightly larger than the plug area as explained in Appendix A. Extrapo- 
lation of the curve to zero heat transported gives an intercept which is 
the difference between an "effective" recovery temperature for a uniform 
temperature body and stagnation temperature. A typical test curve is 
shown in figure 4. 

The test conditions are summarized in table II. In all cases the 
plate leading edge was unswept, and the surface on which the measurements 
were made was at zero angle of attack. 


EXPERIMENTAL RESULTS 
Surface Pressure 


The ratios of measured surface pressure to free-stream static 
pressure, p/p^, are shown in figure 5 as a function of distance x from 
the leading edge. Spanwise surface-pressure variations were found to be 
less than 2 percent within 1/2 inch of midspan. The pressures for plates 
with cylindrical blunting are three to four times larger than free-stream 
pressure near the leading edge, and decrease with distance x. It is to 
be noted that, in contrast, the pressures measured on the square leading- 
edge plate are low near the shoulder, and increase with distance x. 
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However, the pressure variation aft of 1-1/2 inches coincides with the 
distribution of the cylindrical nosed plates. The low pressures are due 
to overexpansion at the square leading edge and have been observed pre- 
viously (refs. 5 and 14) . The present results corroborate the results of 
references 5 and 14 that the overexpansion does not affect the pressure 
distribution aft of about two leading-edge thicknesses. For both types 
of leading edges , the pressures approach free-stream values near the rear 
of the bodies. 


Flow-Field Surveys 


Impact pressures were measured in the flow field for leading-edge 
Reynolds numbers of I 65 O and 6600. These surveys were made to obtain an 
estimate of the boundary- layer thickness, flow quantities along the 
boundary- layer edge, and the shock-wave location. 

In figure 6, the variation of impact pressure with distance above 
the plate is shown for a typical survey. Three distinct flow regions are 
evident in this figure. The first, near the surface, corresponds to the 
boundary layer where the impact pressure increases rapidly. A second 
region is observed to exist just outside the boundary layer where the 
impact pressure is nearly constant. The third region appears to be a 
region of possible shear flow where the impact pressure increases con- 
tinuously with height up to the shock wave. The maximum or peak impact 
pressure corresponding to the shock wave is not shown in figure 6. Above 
the shock wave, the impact pressure drops to the free-stream value. 

The height above the plate at which the peak impact pressure was 
measured was defined as the location of the shock wave. In figure 7> 
the location of the shock wave is plotted versus distance from the lead- 
ing edge. Blast wave theory (refs. 7 and 9) predicts the shock-wave loca- 
tion for regions not too close to the leading edge of blunt flat plates 
in hypersonic flow to be: 

L I (d/ £l = (constant) ( (l) 
d V d / 


This relation, adjusted to fit the data far from the cylindrical leading 
edge, is presented in figure 7 as the solid line. The data agree well 
with the prediction up to within 3 diameters from the leading edge. The 
constant, empirically determined, has a value of 1.47- This constant as 
given in reference 9 is 0.89(C]}) 1/3 . The value calculated from this 
relation for assumed Cj) of 1.2 differs from the experimentally 
determined constant by a factor of 2 2//s . 
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Heat Transfer 


In figure 8, the local "effective" temperature recovery factor is 
shown plotted as a function of distance from the leading edge. The 
recovery factors vary from a value of about O .65 at the shoulder to 
approximately values of 0.8 and 0.86 at x/d = 2, and are nearly constant 
at these values for the remainder of the plate length. It is believed 
that the low values at the shoulder are due to the linear extrapolation 
to zero heat-transfer rate as shown in figure 4 and to the implied sur- 
face condition of uniform body temperature. It is noted that the recov- 
ery factors are higher for the lower Reynolds number conditions. This 
may be experimental error in evaluating the recovery temperature because 
of the corresponding low values of heat-transfer coefficient. 

The measured heat-transfer coefficients are plotted as a function 
of x/d as shown in figure 9- The heat-transfer coefficients are high 
near the leading edge and decrease with distance from the leading edge. 

It is also noted that the heat-transfer coefficients increase with increase 
in free-stream Reynolds number based on leading-edge thickness. 


The edge of the boundary layer was determined from the velocity 
distributions through the boundary layer. The velocities were computed 
from the impact pressures on the assumption that the total temperature 
was constant through the flow field and the static pressure gradient 
normal to the plate was equal to zero. The edge of the boundary layer 
is defined as that distance above the plate surface where the velocity 
ratio has approached within 1 percent of the straight "reference" line 
drawn through the experimental points in region II , see figure 10 . 

The resultant boundary- layer thicknesses are presented in figure 11. 
The theoretical prediction of boundary- layer thickness on a flat plate 
with zero thickness leading edge, reference 1, as altered by introduction 
of the leading-edge thickness, d, is given below: 


DISCUSSION 


Boundary-Layer Thickness 


6 

d 



4.27 ~m» 2 



( 2 ) 


Moo - n/ Recod 


Equation (2) is compared with two sets of data in figure 11 for cylindri- 
cal blunting. It is noted that the boundary- layer thickness for the thin 
leading-edge plate is predicted fairly well by equation (2) , except for 
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the region near the leading edge. The boundary- layer thickness measured 
over the thick leading-edge plate is somewhat lower than predicted by 
the theory. 


Surface Pressure Distribution 


References 5 and 6 indicate that a linear combination of a viscous 
term and an inviscid term may be the proper type of formulation to 
describe the surface pressure distribution over blunt flat plates in 
hypersonic flow. The following linear combination of pressure terms is 
proposed: 


JL = 1 + abX + Pci 


(3) 


where a and 3 are empirical constants to be determined by comparison 
with experiment. 


The viscous pressure term, bX, is the parameter pointed out in ref' 
erence 1 as a result of the theoretical study of hypersonic flow over a 
flat plate with a zero thickness leading edge, where. 


X 



b 


"0.865 


Moo 2 

vw 


0.166(7-1) 


7 


(4a) 

(4b) 


= Pw ^ 

M-oo T w 


(4c) 


The inviscid pressure term, cl, is obtained from the 'blast wave 
analogy to the hypersonic flow over blunt flat plates, references 7 and 8, 
and also from a characteristic development of the flow field described 
briefly in reference 9> where 


I 



(5a) 


The factor c in the inviscid pressure term, as explained in references 
8 and 9, is proportional to the energy released into the transverse flow 
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field by the blunt leading edge. The more exact formulation of this 
factor, reference 9 > is given as 

c = Cl (C D ) 2/3 (5b) 


The value of c x in the above equation is given in reference 9 as 0.112 
for air and O.I 69 for helium. The pressure drag coefficient, Cp, for 
the cylindrical leading edge is approximately 1.2 (see, e.g., ref. 9 ) • 

For a square leading edge the drag coefficient is found to be approximately 
1.8 if one assumes impact pressure acting over the complete front face. 

A correlation of the surface pressures on flat plates in supersonic 
and hypersonic flow is shown in figure 12 where measured pressures are 
compared with values obtained from the following form of equation ( 3 ) : 


2 - = 1 + _b*_ + (6) 

Poo ^ 2 2/3 


It has been found by trial that values of 1/ n/ C w for a and l/ 2 2/3 for (3 
result in the best correlation of the data of the present tests, as well 
as data from references 4, 5* 6 , 15, and 1 6 , see table III. The solid 
line represents a 1 -to-l correspondence of measured to theoretical values. 
The data are noted to be within approximately 20 percent of the solid 
(theoretical) line. The calculation has been made for a temperature 
recovery factor assumed equal to 0.86 for the data of references 4 , 5 , 
and the present tests. The recovery factor for the helium data of refer- 
ence 5 was assumed to be 0.82. Calculations for the conditions of refer- 
ences 15 and 16 were made using their measured wall temperatures. The 
wall temperature affects the parameter b and Cw of the viscous term 
in equation ( 6 ) . In view of the considerable range of variables encom- 
passed by this correlation, the semiempirical expression in equation ( 6 ) 
is believed to be quite general in applicability. 


Local Mach Number 


The local Mach number along the boundary- layer edge was obtained by 
two methods which are compared in figure 13 . Both methods utilize the 
measured surface pressures and the assumption of zero pressure gradient 
normal to the plate. The set of points in figure 13 labeled "local 
measurements" was obtained by use of local measured impact pressures. 

The second set, labeled "reduced total pressure," was obtained by assum- 
ing the total pressure along the boundary- layer edge to be constant and 
equal to the total pressure behind the leading-edge normal shock wave. 
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Two results are noted from figure 13 . First, although the Mach number 
is seen to increase slowly with distance from the nose, it does not 
reach free-stream Mach number within the test length. The second and 
more important result is the excellent agreement between the Mach number 
distributions obtained by the two different methods. 

Reduced local Mach numbers have been reported previously by Brinich 
(ref. 17 ) as obtained in tests conducted at a Mach number of 3.1 over a 
hollow cylinder. Both sharp and square leading edges were tested, and 
the local Mach number was found to be reduced up to a distance of approx- 
imately 3000 nose thicknesses back from the leading edge. Also, Crawford 
and McCauley have reported (ref. 18) reduced local Mach numbers obtained 
over a hemispherical-nosed cylinder at free-stream Mach number of 6.8 for 
distances up to 3 nose thicknesses back along the afterbody. 

From the results of the present tests and those of references 17 
and 18 , it can be concluded that the local properties along the boundary- 
layer edge can be obtained by the reduced total pressure method for the 
range of the present tests and for similar tests of other investigations, 
if a detached normal shock exists at the leading edge. 


The assumption that the total pressure was reduced and constant 
along the boundary-layer edge was applied to the present data, and the 
local Reynolds numbers were computed along the boundary- layer edge. In 
figure Ik the ratio Rex/Re^ is shown plotted as a function of dimen- 
sionless distance, x/d. The quantity Re x is based on local quantities 
at the edge of the boundary layer and distance x. The quantity Re^ 
is based on free-stream conditions ahead of the leading- edge shock wave 
and the diameter of the leading edge. The data are linear on the loga- 
rithmic plot. The best fit to the data, determined by method of least 
squares, resulted in the relation 


This empirical relation is shown as the solid line in figure 14. It is 
noted that the local Reynolds number at the boundary- layer edge is 
decreased to about 40 percent of the free-stream value and this percent- 
age decrease does not vary for changes of 10 to 1 in leading-edge Reynolds 
numbers. Reduction in local Reynolds number has also been obtained by 
Brinich (ref. 17 ) in his transition studies. 


Local Reynolds Number 



(7) 
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In view of the successful correlation of the present data in terms 
of local Reynolds number, the method described above was applied to the 
results of Kendall (ref. 4), Hammitt, Vas, and Bogdonoff (ref. 5), Bertram 
(refs. 6 and 15), and Erickson (ref. 16) . The results are presented in 
figure 15. The solid line represents a datum line or variation where the 
local Reynolds number per inch is equal to free-stream Reynolds number 
per inch. It may be noted that in figure 15 the effects of leading-edge 
thickness and free-stream Reynolds number, for each set of points, are 
fully accounted for. Also in figure 15 the Reynolds number ratio for a 
fixed value of x/d decreases with increasing Mach number and seems to 
approach a limiting value as Mach number increases. There are two other 
parameters which differ among these sets of data, namely, leading-edge 
shape and test gas. In the previous section on surface pressures, the 
shape of the leading edge has been shown not to affect the surface pres- 
sure distribution past an x/d of 1 to 2. Thus we may rule out the 
leading-edge shape parameter provided a detached shock exists, leaving 
Mach number and ratio of specific heats as possible parameters affecting 
correlation of the sets of data. 

The empirical fit to each set of data in figure 15 was found by. the 
method of least squares. The general relationship is represented by the 
following equation. 


The values of K and a are tabulated in table IV together with those 
for the present tests. 

The slope is practically the same for all the sets of data shown 
in figure 15. Thus the test gas and Mach number dependence is contained 
in the parameter K of equation (8) ■ An analysis presented in Appendix B 
suggests a correlation parameter, namely, the total pressure ratio across 
the normal shock occurring at free-stream Mach number. Thus the values 
of the parameter K were plotted versus the total pressure ratio across 
the normal shock as shown in figure 16. It is noted that the parameter K 
is correlated well by the following equation. 


An average value of the slope, a, was taken to be 0.91. This average 
value of a and the relation for K in equation (9) when substituted 
in equation (8) give the best fit to all the data sets. 



( 8 ) 



(9) 
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( 10 ) 


The excellent correlation obtained hy the use of equation (10) is 
shown in figure 17. It should be noted that the data correlated by use 
of equation (10) cover an extremely wide range of variables. For exam- 
ple, the data include a Mach number range of from 4 to 17.3, leading-edge 
thicknesses of 0.0002 to 1.0 inch, and free-stream Reynolds numbers per 
inch from 2500 to 10 6 . In addition, the data include tests both in air 
and in helium. The success of this correlation suggests that in all cases 
considered the total pressure along the boundary- layer edge is constant 
and equal to the total pressure behind the leading-edge normal shock. 
Further tests are required to determine the applicability of the correla- 
tion over ranges of variables beyond those covered by the data shown. 


Heat Transfer 


The heat-transfer coefficients are presented in dimensionless form 
in figure 18. where the local Nusselt number is plotted as a function of 
the local Reynolds number. The well-known Pohlhausen solution for local 
heat transfer in the laminar boundary layer over a flat plate is 


Nu x = 0.332(Pr) l/3 


(ID 


Because the Prandtl number of the present tests is nearly 0.72, equa- 
tion (ll) reduces to 


Nu x = 0.295 \/Re x 


( 12 ) 


It is to he noted in figure 18 that this simple Pohlhausen solution fits 
the majority of the data quite well. The scatter from the line of equa- 
tion (12) can he explained hy the limitations of the present heat-transfer 
instrumentation. The heat-transfer elements yield inaccurate results 
when h falls to values of around 0.5 Btu/ft 2 , hr, °F. The tare then 
Becomes equal or greater than the net heat transfer to the stream. 

The effect of blunting on the heat-transfer coefficients over a flat 
plate can now he examined. The correlation of figure 18 was made on the 
basis of properties evaluated at local conditions at the edge of the 
boundary layer. However, in order to account for nonlinear variation 
of viscosity and thermal conductivity, the well-known T’ method 
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(ref. 19) for evaluating properties was incorporated. As shown in 
Appendix C, the ratio of heat-transfer coefficients for blunt and sharp 
plates is 



Note that the ratio of Reynolds numbers appearing under the square root 
in equation ( 13 ) can be obtained directly from the present correlation 
as given by equation (10) . The primed quantities are evaluated from the 
equations of Appendixes B and C and the known pressure distribution over 
the surface of the blunted plate. The quantities applying to the sharp 
plate are evaluated from free-stream conditions. 

The results are now applied to a flat plate in flight for the fol- 
lowing set of assumed conditions: 


K, 

= 10 

Re OT /ft 

= 10 6 

Too 

= 420° R 

Tv 

= 1000° R 

d 

a 1 in. 


The reduced total pressure is assumed to persist over the entire 
surface. The local static pressure and local Mach number distributions 
are presented in figures 19 and 20 for these conditions. The pressure 
was calculated by equation ( 3 ), and the local Mach number was calculated 
from equation (B5) . The ratio of the heat-transfer coefficients for a 
blunted plate to the heat-transfer coefficients for a sharp plate was 
calculated by equation (13). The variation of this calculated ratio with 
distance from the leading edge is shown in figure 21 . 

It is noted from figures 19, 20, and 21 that the effect of blunting 
the leading edge is to increase the heat-transfer coefficients in the 
region over the plate surface where high static pressures exist. How- 
ever, as noted from figure 21 , the effect of blunting the leading edge 
is to reduce the heat-transfer coefficient far back on the plate, where 
surface pressures have a value near that of the free stream. It is thus 
concluded that blunting with consequent reduction of total pressure along 
the boundary- layer edge increases the heat-transfer coefficient near the 
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leading edge where the high static pressures exist, and reduces the 
heat-transfer coefficient far hack on the plate where the static surface 
pressure approaches the free-stream value. 


SUMMARY OF RESULTS 


Pressures and heat-transfer rates -were measured on hlunt flat plates 
at a nominal Mach number of 4 for free-stream Reynolds numbers per inch 
of 2380 and 6600 in air. The Reynolds number based on leading-edge thick- 
ness ranged from 600 to 6600 . Surface pressures near the leading edge 
were found to be nearly three to four times larger than free-stream static 
pressures, and these higher pressures persisted far back on the plates. 

A linear combination of the viscous hypersonic parameter and the blast 
wave inviscid parameter was found to correlate pressure distributions 
obtained in this and other investigations for Mach numbers of 4 to 17 . 3 * 
free-stream Reynolds numbers per inch from 2500 to 10 6 , and leading-edge 
thicknesses from 0.0002 to 1.0 inch. 

Boundary- layer thicknesses, obtained for leading-edge Reynolds num- 
ber of I 65 O and 6600, agreed with the theory for flat plates with zero 
thickness leading edge. Experimentally determined flow quantities along 
the boundary- layer edge agreed remarkably well with values calculated 
with the assumption of reduced total pressure equal to the leading-edge 
stagnation pressure. 

The local Reynolds number at the boundary- layer edge was found to 
be decreased to about 40 percent of the free-stream value by blunting 
at a Mach number of 4. The decrease in local Reynolds number was not 
appreciably affected by changes of 10 to 1 in leading-edge Reynolds 
number.' 

A general correlation of the ratio of local Reynolds number to 
leading-edge Reynolds number was obtained by comparing the present test 
results with the results of similar tests which covered wide ranges of 
Mach number, leading-edge thickness, free-stream Reynolds number, and 
test gas. The basic assumption under lying this correlation is that the 
total pressure along the boundary- layer edge does not differ appreciably 
from stagnation point pressure. If this condition is assumed to be true, 
then the local Reynolds number along the boundary- layer edge is lower 
than free-stream Reynolds number by a factor which is very nearly the 
square root of the total pressure ratio across the normal shock wave at 
the leading edge. 

Heat-transfer rates were measured for leading-edge Reynolds numbers 
of 238 O and 6600. Local Nusselt numbers were correlated by local Reynolds 
numbers and predicted by the well-known Pohlhausen solution wherein local 
conditions at the boundary- layer edge are used. As compared to the sharp 
condition, blunting of a flat plate was found to increase the heat-transfer 
coefficients in the region near the leading edge where high static 
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pressures and reduced total pressure exist, and to decrease the heat- 
transfer coefficients in the region far from the leading edge where 
reduced static and total pressures exist. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 9, 1957 
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APPENDIX A 

DETERMINATION OF THE HEAT -TRANSFER AREA 


The net heat was considered to be transferred from the top surface 
of the cylindrical test plug to the air stream. Thus the area, S, used 
in the following equation was the area of the top of the test plug. 


q = hS(T-T n ) 


(Al) 


However, the film stretched over the model (see fig. 3 insert) does con- 
duct some heat away from the plug. Also, the film receives energy from 
the cylindrical surface of the plug by free molecular conduction through 
the annular air space around the plug. This film acts much as a circular 
fin in dissipating the heat from the test plug. This effect is to 
increase the area to be used in calculating the heat-transfer coefficient. 
The following analysis of the fin 
to the test area. 

The differential equation 
governing the temperature distri- 
bution in the circular film fin 
may be found by summing the quan- 
tities of heat transferred by the 
various means to and from an 
annular element of the fin (see 
sketch (a)). Azimuthal varia- 
tions of these quantities around 
the test plug will be considered 
negligible. The assumption is 
made that temperature differences 
are small, so that the radiation 
exchange terms may be written in 
linear form. The width of the 
annular air space is of the order 
of a mean free path of the gas. 

Thus, it is assumed that the 
annular element of. fin gains heat 
from the plug by free molecular conduction through the annular air space. 
This element also transfers heat to the stream by convection through the 
flow boundary layer. A further assumption will be made that the varia- 
tion of heat-transfer coefficient and recovery temperature is negligible 
over the area of film considered. When these terms are combined with 
the normal terms of conduction in the film fin, the heat flow balance 
may be written for the annular element 

dqo + dqc + dqp - dq, ± - dq^ = 0 


effect leads to a form of the correction 



d 9fm 

Sketch (a) 
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The differential equation is then found to he: 


dr" 


T - - ) + i — (t--)-b(t--) = 0 
B/ r dr \ By V B 


(A2) 


The differential equation (A2) is a form of Bessel* s equation. The 
solution may he written in terms of modified Bessel functions of zero 
order , first and second kinds , as 


T - - = AgloCrN/B) + A 4 Ko(rN/B) 
B 


(A3) 


where 


A = 


_ ( 3M m K 

2 H. 


t-pk 


fKf 


+ 40eT o 3 — + h 

Tr, T-n 


(Ah) 


B = 


3NV mK 


tfkf \ 2 \frt 


+ 4aeT 0 3 + h 


(A5) 


and Aq and A 4 are constants of integration to he determined hy the 
following boundary conditions at 


r = r p , T = T p 


r = r w 


T = T w = T- 


(AS) 
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In the range of Interest, the modified Bessel functions in the 
solution (A 3 ) may be replaced by the asymptotic expressions, reference 20, 
for large values of the argument (r'/l). These expressions are 


I 0 (r n/b) = 

v2n(r n/b) 


K 0 (r n/b) = / 5 exp(-r *>/B) 

v 2 (r n/b) 




(A7) 




Relations (A 7 ) are introduced into equation (A 3 ) , the boundary conditions 
are applied, and the constants of integration are evaluated. The radius, 
r m , at which the minimum temperature occurs is found to he very nearly 
the average radius given by 


r m 


*w + r p 

2 


(A8) 


The temperature distribution in the circular fin is then given by the 
resulting form of equation (A3) as 




exp(r VB) + ^ 

2r m '/B-l\ 
2 r m -s/Bfl/ 

exp( 2 r m s/B-r n/b) 

exp(r p n/b) + 

/ 2r m n/B-1' 
^ 2 r m n/b+ 1 ' 

) ex P ( 2r m '/B-r p '/!) 


(A9) 


The heat removed from the plug by the film is transferred to the stream 
according to the following relation 


dq = h(T-T r) )dS 


(A10) 


This expression is integrated using the temperature distribution in 
equation (A9) for T, and assuming a constant heat-transfer coefficient. 
The resulting expression is found 
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<1 


p T m 

2ith J (T-T n )r dr 


(All). 


Assuming that the amount of heat represented in equation (All) were to 
be transferred at plug temperature from an area given by kr 2 , we have 

q = rth^-rp 2 ) (Tp-T^) (A12) 


Thus by equating the right-hand sides of equations (All) and (A12) we 
can express the heat-transfer radius r as follows: 



In actual computation of the correction, the emissivity of the film 
was assumed to be the same as that of the chrome-plated plug surface, 
that is 0.1. The emissivity of the plug with film was found to check 
closely with the value normally taken for polished chrome. The conduc- 
tivity of the film was taken as 0.1 Btu/hr, ft 2 , °F/ft. This value was 
obtained from manufacturer’s literature, and was not checked experimen- 
tally during these tests. The first approximation to h, found by 
using rp, was used to determine r. The correct heat-transfer 
coefficient is then found using this r computed from equation (A13) . 
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APPENDIX B 

LOCAL REYNOLDS NUMBERS 


The ratio of local to free-stream Reynolds number is analyzed for 
the case of a blunt flat plate in supersonic flow. The general assump- 
tion is made that the total pressure along the boundary- layer edge is 
lower than the total pressure in the free stream ahead of the bow shock 
wave. In addition, the static pressure is assumed constant through the 
boundary layer. The local Mach number at the boundary- layer edge is 
specified by the local static and total pressures. The local static 
temperature can be found from the local Mach number and the known 
stagnation conditions. 

The defining expression for the ratio of Reynolds numbers is given 
below . 


Rex 

Re<xd 



(Bl) 


In order to relate the various terms in equation (Bl) to the ratio 
of local static to total pressure, the following equations are obtained 
from reference 21. 




(B2) 





(B3) 


The viscosity is assumed to vary as a power of the absolute temperature 


H~(T) W 


(B4) 


The following ratios are easily obtained from equations (B2), (B3) , 
and (Bk) . 
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Hx, 


l-(TAt 5 ) 

l-(Tco/T t5 ) 


1/2 


y- 1 

Pti/ 


7-1 

= 

T Vp p tl / 




Hi M 

!^y 

Pti/ . 


(B5) 


(b6) 


(B7) 


The ratios as determined in equations (B5) , (B6) , and (B7) are introduced 
into equation (Bl) . The resulting expression is 


Rex 

Re ood 



1-(T/T t6 ) 

l-(Too/T t5 ) 


1/2 



(B8) 


where the constant 


£ is as follows: 


i = 



(B9) 


If the assumption is made that the local total pressure is reduced 
to that value of total pressure existing behind the leading-edge normal 
shock, the ratio of Reynolds numbers becomes. 


R e x 

Reood 




1-(T/T tl ) 


l-doo/rti) 


1/2 



(BIO) 


The total pressure ratio term is a function only of free-stream Mach 
number and gas. As is shown in the correlations of figure 15, the braced 
term is constant for fixed free-stream Mach number. 
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Prohahle values of the exponent £ may he calculated. For a tem- 
perature range from 3 OO 0 R to 900° R in air, w is nearly 0.75* From 
equation (B9) the value of £ is found to he 0 . 5 . For helium, w is 
approximately O .63 for a temperature range of 10 ° to 100 ° R, and £ is 
found to he O. 65 . 
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APPENDIX C 


REFERENCE TEMPERATURE FOR FLUID PROPERTIES 


The temperature at which fluid properties are evaluated is arbitrary 
in empirical correlations. In reference 19 an expression is obtained 
which gives an evaluation temperature, the well-known T* temperature. 

The advantage of the use of the T* temperature is that when properties 
are evaluated in this manner, the drag and heat transfer can be expressed 
in such a form as to minimize their dependence on Mach number, wall 
temperature, Prandtl number, and power law exponent for viscosity and 
thermal conductivity. The expression derived in reference 19 is presented 
below . 


= 1 + 0.032 M 2 + 0.58 
T 



(01) 


If the Nusselt number, as given by the normal Pohlhausen solution, 
is written in terms of the T 1 properties, the form can be shown to be. 



It may be noted from equation (C2) that the T* method of obtaining 
reference temperature .merely takes into account the nonlinearity of the 
viscosity and thermal conductivity variation with temperature. Thus, 
if an experiment is conducted over a range of temperatures where the 
viscosity and conductivity vary linearly with temperature, results 
obtained for heat-transfer correlation will be unaltered by the T’ 
method of data reduction. 

For the case of blunt leading edges, the properties are evaluated 
at local conditions, the T* method is applied to the local temperatures, 
and equation (C2) becomes: 



For the case of the sharp leading edge, the local properties are 
normally considered as based on free-stream conditions ahead of all dis- 
turbances, and the T* method is applied as follows: 
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h 


t 

s 



0.295 s/iei 


00X 


(ck) 


The following ratio of heat-transfer coefficient is obtained from 
equations (C3) and (C4). 



(C5) 
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TABLE I.- STREAM CONDITIONS 


Static pressure. 

Mach 

Stream diameter. 

microns Hg abs 

number 

in. 

300 

3-95 

3.6 

100 

3-87 

3-0 


TABLE II.- SUMMARY OF PRESENT TEST CONDITIONS 


Type of test 


Model leading edge 


Thickness, 
in. 


Shape 


Test conditions 


Moo 


Reoo/in. 


Figures showing 
data 


Surface pressure 


Local flow field 

Heat-transfer and 
recovery factor 


iA 

1/2 


1 



1 


iA 

i 

I 


Cylindrical 

I 

Square 
Cylindrical 

1 

Square 
Cylindrical 


v 


3-95 


3-87 


Y 

3-95 

3-95 

3-87 


6600 


Y 

2380 


v 

6600 

I 

66oo 

2380 


5,12,14,15,17 


7,H,13 
6,7,10,11,13 
8,9,18 


i 


TABLE III.- TEST CONDITIONS OF OTHER INVESTIGATIONS 


Investigation 

Re ood 

Mco 

Gas 

Leading edge 

Kendall (ref. 4) 

12 to 50 

5*8 

Air 

Sharp 

Bertram (ref. 6) 

370 to 1900 

6.86 

Air 

Square 

Bertram (ref. 15) 

38 to 55 

9.6 

Air 

Sharp 

Hammitt and Bogdonoff (ref. 3) 

3480 to 6100 

11.8 

Helium 

Square 

Hammitt, Vas, and Bogdonoff (ref. 5 ) 

15,000 to 38,000 

12.7 

Helium 

Square 

Erickson (ref. 1 6 ) 

650 and 685 

16 and 17.3 

Helium 

Sharp 


TABLE IV.- REYNOLDS NUMBER CORRELATION PARAMETERS 


Investigation 

Moo 

7 

K 

a 

Present tests 

3-95 

1.4 

0.415 

0.91 

Kendall (ref. 4) 

5.8 

1.4 

.223 

.95 

Bertram (ref. 6) 

6.86 

1.4 

.152 

.93 

Bertram (ref. 15 ) 

9.6 

1.4 

.117 

.87 

Hammitt and Bogdonoff (ref. 3 ) 

11.8 

1.67 

.129 

.89 

Hammitt, Vas, and Bogdonoff (ref. 5 ) 

12.7 

1.67 

.132 

.88 

Erickson (ref. 16) 

16 

1.67 

.096 

.90 

Erickson (ref. 16) 

17.3 

1.67 

.076 

.93 
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Figure 1.- General arrangement of wind-tunnel test section. 
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Figure 2.- Variation of Mach number with distance from center of stream 
for Reynolds number per inch of 66 00 at axial distance of 1.25 inches 
from exit plane of nozzle. 
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Figure 3.- Heat-transfer tody 
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Figure 4.- Variation of heat -transfer rate with difference between sur- 
face temperature and stream stagnation temperature for free-stream 
conditions of Moo = 3-95 and Re^ = 66 00. 
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Figure Variation of the ratio of surface pressure to free- stream 
pressure with distance from the leading edge. 
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Figure 6.- Variation of the ratio of impact pressure to leading-edge 
stagnation pressure with distance from the plate surface. 
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Figure 7.- Variation of shock-wave location with distance from leading 
edge for = 3.95 and Re^ per inch of 6600. 
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Figure 8.- Variation of the local "effective" temperature recovery factor with distance from the 

leading edge. 
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Figure 9 .- Variation of heat-transfer coefficients with distance from 
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Figure 10.- Variation of the velocity ratio through the boundary layer. 
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Figure 11.- Comparison of measured boundary- layer thickness with theory. 
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Comparison of measured pressure ratios with calculated pressure ratios 
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( a ) Re^ = 6600 

Figure 13 .- Variation of Mach number along the ho undary- layer edge. 
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(b) Re^ = 1650 
Figure 13 •- Concluded. 
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Figure 14.- Variation of the ratio of local Reynolds number at the boundary- layer edge to the 
free- stream Reynolds number with distance from the leading edge. 





NACA TN klb2 





Reference 


" Symbol Mg. 

3.87 

|s o and Present tests 

1 3.95 

□ 5.8 Kendall 

O 6.86 Bertram 

V 96 Bertram 

| A 118 *HVB 

| t, 127 HVB 

lb 16.0 Erickson 
A Q 17.3 Erickson 

tTHammitt, Vas, and Bogdonoff 


Distance from leading edge, x/d 

Figure r?.- Comparison of Reynolds number ratio from various sources, see table III. 
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Figure l6.- Variation of the parameter K from various sources with total pressure ratio across 

a normal leading-edge shock wave. 
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Figure 17 . - Correlation of results calculated from present tests and 

from other sources. 
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Figure 18.- Variation of local Nusselt number with local Reynolds number for a cylindrically 

blunted flat plate in supersonic rarefied air flow. 
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Figure 19. - Surface pressures calculated for a cylindrically blunted flat wing for zero angles 

of yaw and attack at an altitude of 100,000 feet. 
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Figure 20.- Mach number ratio calculated for a cylindrically blunted flat wing for zero angles 

of attack and yaw at an altitude of 100,000 feet. 
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Figure 21.- Effect of "blunting on the heat-transfer coefficients 
zero angles of attack and yaw at an altitude of 


calculated for a flat wing for 
100,000 feet. 
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